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Tomography, in most formulations, 
requires an incoherent signal. For a con- 
ventional transmission electron micro- 
scope, the coherence of the beam often 
results in diffraction effects that limit the 
ability to perform a 3D reconstruction 
from a tilt series with conventional tomo- 
graphic reconstruction algorithms. In this 
paper, an analytic solution is given to a 
scanned Gaussian beam, which reduces the 
beam coherence to be effectively incoher- 
ent for medium-size (of order 100 voxels 
thick) tomographic applications. The 
scanned Gaussian beam leads to more 
incoherence than hollow-cone illumina- 
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1. Introduction 

To perform microtomography, the first question is 
whether to use a coherent technique such as holoto- 
mography [1], or an incoherent technique, which ties 
the signals to mainstream tomography [2]. Amorphous 
samples are usually regarded as having incoherent 
interactions even with coherent beams [3], although 
coherent techniques have been applied to amorphous 
samples [4]. In biology, the samples are intrinsically 
amorphous if the growth of ice crystals is avoided. The 
now-standard technique of plunge freezing inhibits the 
formation of ice crystals and allows the preparation of 
amorphous samples. Tremendous progress has been 
made [5]. 

If the sample is crystalline or contains a significant 
proportion of crystallites, the incoherence must come 
from elsewhere. A suitable detector may observe an 
incoherent signal, as in the High-Angular Aperture 



Dark Field (HAADF) method, which exploits the fact 
that random motion of the nuclei in a crystal lead to an 
incoherent signal at large scattering angles [6]. (In the 
theoretical description, the Debye-Waller factors 
become very small.) Recently, the complement of the 
ideal HAADF signal has been used, under the name of 
Incoherent Bright Field (IBF) imaging [7] with the 
additional advantage that the signal is a monotonic 
function of thickness whereas a practical HAADF 
detector (with an inner diameter and an outer diameter) 
has a single-peaked thickness-intensity relation. 

In this paper, a different approach is taken. A propos- 
al is made to scan the beam in a particular manner 
while simultaneously recording an image to create a 
signal that is effectively incoherent. I call this approach 
"synthetic incoherence" to distinguish it from the more 
common situation of incoherence induced by thermal 
fluctuations of the source [8]. Using the definition of 
coherence, it is shown that it is possible to reduce the 
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coherence of a beam across a surface to a level which 
is negligible tomographically. This is accomplished by 
scanning an out-of-focus Gaussian beam across the 
sample using a position-varying dwell time that itself 
has a Gaussian distribution. Of course, an incoherent 
beam cannot be strictly parallel, and this lack of paral- 
lelism sets a limit for the thickness of the sample 
assuming that the projection assumption [2] is required 
for the reconstruction, i.e., assuming that the beam trav- 
els in a straight line. 

The definition of the mutual coherence function is 
given in many sources, such as Spence [9]; 

T l2 (r l ,r 2 ,T) = \\m^-[ dt\j/*{f l ,i)Xj/{f 2 ,t + T) (1) 



where r x and r 2 are two positions in the optical field y/, 
t is a time, and T is a time difference. Only the equal 
time case, i.e., T= 0, is of interest here [10]. The com- 
plex degree of coherence is defined as the normalized 
mutual coherence function 



Yn(n,r 2 )- 



T u(fi,r 2 ) 



[r„©r 22 (F 2 )] 1/2 



(2) 



The normalized mutual coherence function is an inten- 
sity-weighted average of the phase factor between two 
illuminated points. When |y 12 | « 1 the beam is incoher- 
ent. Of course, y 12 (r, r) = 1, so the best one can do is to 
have \Yi 2 (?i, r 2 )\ go to zero rapidly as |?i - f 2 \ grows. 



2. Gaussian-Scanned Gaussian Beams 

The Gaussian wave function solution to the paraxial 
ray approximation for a beam traveling with its center 
in the z direction is [11] 



y/(x,y,z)- 



1 



1 + iZ 



exp 



X 2 + Y 2 
1 + iZ 



(3) 



where the scaled dimensions X, Y, and Z are related to 
the physical dimensions x, y, and z by 



(4) 
(5) 
(6) 



Z = 



'o 
Xz 
nr 2 



where r is a parameter characteristic of the beam waist 
corresponding to two standard deviations of the inten- 
sity, i.e., the square of the wave function and X is the 
wavelength. The beam waist may be given in terms of 
the angle Q y between the direction of maximum intensi- 
ty and the angle at which the intensity has fallen by a 
factor of e 2 



/rft 



(V) 



The time average of Eq. (2) may be performed by vary- 
ing the beam conditions in a controlled fashion. This is 
in contrast to the usual practice of performing the time 
average of a source subject to statistical fluctuations. 
An analytic solution is achieved in the following case. 
The beam is scanned in the focal plane defined by z = 
0. The beam center is given by (x ,y ). We restrict the 
discussion to the plane of a surface of the sample 
defined by the coordinate z, i.e., a plane whose normal 
is parallel to the beam direction. The mutual coherence 
function in the plane is given by 



T i2 (x l ,y i ,z,x 2 ,y 2 ,z)=\ J dX^dY — 



\+Z L 



exp 



(X-Xyf+jY-Yj (X 2 -X f + (Y-Yj 



1-iZ 
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<X ,Y ) 
(8) 



where w(X , Y ) is a weighing factor which may be 
taken as proportional to the dwell time. The normaliza- 
tion 



1: 



- jldX dY w(X ,Y ) 



(9) 



is imposed. Assume that w(X , Y ) = w(X )w(Y ), i.e., w 
is a function of two variables factorized into two func- 
tions of one variable (which happen to be the same). 
The one-dimensional normalization 



1: 



:jldX W(X ) 



(10) 



is imposed, which implies the two-dimensional normal- 
ization of Eq. (9). 
The function 



f{X„X 2 ) 
= \ dX exp 



/ 



(x.-x,) 2 (x 2 -x )M 



\-iZ \ + iZ 

has an integrand that may be expressed as 



w(X ) (11) 
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/; = exp 



X\ + X\ - 2(X 1 + X 2 )X +2Xl ^ 



1 + Z 2 



exp| ,"— —[-X 2 +X 2 2 +2(X 1 -X 2 )X ] \w(X ). (12) 



The factor not dependent on X is 

Z 



C = exp 



x 2 +x 2 ^ 
i+z 2 



exp i 



\ + Z' 



(-X 2 +X 2 ) (13) 



Specialize the weight to a Gaussian 



w(X ) = 



2n I 



exp 



( xC 

2V 



(14) 



where X = olr is the scaled standard deviation and cris 
the dimensional standard deviation. For this weight 
function 

f(X„X 2 ) = -^£ dX cxp(^xX 2 ) 



exp 



2k 1' 

-2(X 1 +X 2 )+2iZ(X 1 -X 2 ) v 



(15) 



where a = 2(1 + Z 2 )" 1 + (2Z 2 )" 1 . This integral may be 
performed analytically using the result 

i dxexp(-ax 2 +bx) = J— exp — . (16) 

J — va I 4a I 



So 



f(X 1 ,X 2 ) = 
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exp 
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The ID normalization constant may be obtained by set- 
ting X 2 to X x in Eq. (17), then returning to Eq. (17) and 
setting X x \oX 2 , multiplying, and taking the square root. 
The result is 
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The quotient simplifies to 
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exp 
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Note that J n {f y , r 2 ) at z : =z 2 = z is simply Eq. (19) 
times Eq. (19) under the substitution (X h X 2 ) -> (7 b 7 2 ). 
Let R 2 = {X, -X 2 f + (7 : - Y 2 f. Finally, 



( 



Yn 1= ex P 



R 2 
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= exp 



2I 2 R 2 
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Note that |/ 12 | depends on the coordinates only through 
the difference variable R. Eq. (20) obeys y V2 =1 for R = 
0, as required. If £ — > °°, |7 12 | = exp(-i? 2 /2). This means 
that if the beam is scanned uniformly over the entire 
(infinite) focal plane, the range of the mutual coherence 
is limited to about the focal spot size. (R is the distance 
between any two points in the plane in units of the focal 
spot size parameter r .) On the other hand, if there is no 
scanning, i.e., if £ — > 0, then \y n \ = 1 confirming that a 
single Gaussian is perfectly coherent. 

How narrow can the mutual coherence function be 
made in practice? The intermediate cases may be para- 
meterized by a constant c x defined by c\ = (1 + Z 2 )/Z 2 , 
in which case the coherence scale length is sj\ + c 2 1 4 
times the asymptotic limit. The case c 2 = 12, for exam- 
ple, is a reasonable compromise which achieves a mod- 
erate scan range and for which the scale length of the 
normalized mutual coherence function is twice its 
asymptotic minimum. To simplify the analysis, assume 
Z» 1, which will be true for typical TEM illumina- 
tion. Then we may approximate CjS = Z or 



a = 



zft 



(21) 



in the unsealed variables. Now z6 { is the characteristic 
scale length of the illuminated region, which is about 
the sample size. Because we have a Gaussian weighting 
function, we make a very small error by omitting the 
tails. In practice, the scanning must take place over a 
finite range, which we may call c 2 a. The case c 2 = 6 
represents the omission of the tails past 3 standard devi- 
ations; the total scanned length is V3z0j in this case, i.e., 
roughly twice the sample size, which is a practical 
result. 

The situation imposes a limit for projection tomogra- 
phy. A typical figure for 6 { is 3 mrad. Hence, the illumi- 
nation strikes the sample by rays which are going in 
different directions on this order. If the results are to be 
applied in projection tomography, the number of pixels 
must be below some constant or order 1 over 6 { , which 
means below some constant times 333 for the case of 
#! = 3 mrad. In practice, tomographic reconstructions 
are performed over 100 to 1000 voxels, so the require- 
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ment is consistent with smaller reconstructions but not 
larger ones. 

3. Hollow-Cone Illumination 

Hollow-cone illumination is frequently used to 
enhance the incoherence of the beam [12]. Several 
electron microscope vendors allow such a mode in 
hardware. The analysis of hollow-cone illumination is 
deferred to a future paper [13]; here, a numerical exam- 
ple is provided. In hollow-cone illumination with a 
Gaussian beam, the axis of the fundamental solution 
given in Eq. (3) is tipped by the polar angle and 
azimuthal angle </>. The polar angle is held constant 
and the azimuthal angle is varied over its full range of 

to 2k. An example of the resulting mutual coherence 
function is given in Fig. 1 for a 300 keV electron with 
0=3 mrad, hence a wavelength X = 1 .969 pm, and z = 

1 mm. In the example, the focus of the beam remains at 
the same point in space as the azimuthal angle is var- 
ied. The situation of this example is characteristic: the 
loss of confinement appears to have an envelope of 



l/R — a far slower decay than a Gaussian beam scanned 
with a Gaussian intensity as discussed in the previous 
section and also as shown in the figure. 



4. Conclusions 

To achieve a practically incoherent beam for use in 
electron tomography, a novel illumination scheme is 
proposed that combines elements of STEM and TEM: 
the sample is to be completely illuminated throughout 
the observation while the beam is scanned. Moreover, 
the scanning itself is more concentrated near the center 
so that the intensity is given by a Gaussian with a stan- 
dard deviation under user control. The new scheme is 
compared to hollow-cone illumination and holds prom- 
ise for removing the residual long-range correlations in 
the latter case. The new scheme leads to no substantial 
correlations for separations of 1 nm or more in an 
example chosen with realistic parameters. 

The practical implementation of the scheme requires 
the experimentalists to be able to program the controls 
of a STEM-TEM. Such a group would probably 
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Fig. 1. The mutual coherence function vs. the distance between the two points in the plane R = 
[(x 2 -x l ) + (y 2 - y{) ] lr for the Gaussian beam scanned with a Gaussian intensity profile (dot- 
ted line) and hollow-cone illumination (solid line) with c l = 12, 9 = 3 mrad, A = 1 .969 pm, z = 1 
mm. These parameters yield r = 209 pm. The parameters are defined in the text. 
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involve a collaboration with a vendor or the use of one 
of the limited number of electron microscopes dedicat- 
ed to research into the techniques of microscopy. 
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